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A neutrosophic statistic is a random variable and it has a neutrosophic probability distribution. So, 
in this paper, we introduce the new neutrosophic Birnbaum-Saunders distribution. Some statistical 
properties are derived, using Mathematica 13.1.1 and R-Studio Software. Two different estimation 
methods for parameters estimation are introduced for new distribution: maximum likelihood 
estimation method and Bayesian estimation method. A Monte-Carlo simulation study is used to 
investigate the behavior of parameters estimates of new distribution, compare the performance of 
different estimates, and compare between our distribution and the classical version of Birnbaum- 
Saunders. Finally, study the validity of our new distribution in real life. 
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Normal distribution is considered the most distribution used in our real life. Many new distributions are derived 
from normal distribution using different transformations. Two parameter Birnbaum-Saunders (BS) distribu- 
tion is considered one of these distributions. ‘introduced the BS distribution as a statistical model for fatigue life 
of structures under cyclic stress. In the recent years the BS distribution is used in many fields to its theoretical 
arguments associated with cumulative damage processes, its properties, and its relationship with the normal 
distribution. BS distribution is unimodal, positively skewed also it investigated for applications in engineering 
by many authors see». Also, BS distribution has many applications in other fields such as business, environment 
and medicine see *'°. Also, the BS distribution can be obtained as an approximation of inverse Gaussian (IG) 
distribution see”, it can see equal mixture of an inverse Gaussian and its reciprocal see*’, Many statistical prop- 
erties of BS distribution is studied by many authors such that probability density function pdf, hazard function 
(hf) because it plays an important role in lifetime data see”. 


Definition 1: A random variable X is said to be Birnbaum-Saunders distribution with shape parameter a > 0, 
scale parameter 6 > Oand denoted by X ~ BS(q, f), if the probability density function (pdf) and the cumulative 
distribution (cdf) of X are defined as follows respectively. 


1 -l/x 8B x 3/2+B) 
fas) = a exp (5 (F +2 2) J a x>0, (1) 
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where, Erf is the error function. 
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Figure 1 shows the pdf of BS(q, f) for different values of shape parameter a we can the BS distribution is 
unimodal distribution. Also, Fig. 2 shows the cdf of BS(a, 8) for different values shape parameter aw. Two Figures 
show the changes in the distribution curve when the shape parameter takes different values. 

The hazard function of BS distribution is defined as, 


_ Fass B) 
1 — Fps(x; a, B)’ 
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Figure 3 shows the hf of BS(q, f) for different values of shape parameter w. The Figure shows the changes in 
the distribution curve when the shape parameter takes different values. 

The main reason of choosing BS distribution as a fatigue failure life distribution by ' is known that for the 
analysis fatigue data used any two-parameter distribution such as Weibull, log-normal and gamma distribu- 
tions. For importance of BS distribution, we proposed in this paper a new distribution called neutrosophic 
Birnbaum-Saunders distribution and denoted by NBS(ay, By ).NBS(q, B). In the literature of neutrosophic 
statistics, the start work in neutrosophic statistics is introduced by *° when he showed that the neutrosophic 
logic is more efficient than fuzzy logic. Also, Smarandache”’ present the neutrosophic statistics and showed also 
it is more efficient than classical statistics. Neuterosophic statistics is considered as the generalized of classical 
statistics and it is reduced to classical statistics when imprecise observations in the data. For the efficient of 
neutrosophic statistics see also***°. Many authors introduced the neutrosophic probability distributions such 
as Poisson, exponential, binomial, normal, uniform, Weibull and so on see?”*!~* introduced the neutrosophic 
queueing theory in stochastic modeling.***” and** investigated the neutrosophic time series. Recently, many 
authors studied the neutrosophic random variables see*”“° inserted the new notions on neutrosophic random 
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Figure 1. The pdf of BS(a, 8) for different values of shape parameter a. 
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Figure 2. The cdf of BS(q, 8) for different values of shape parameter a. 
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Figure 3. The hf of BS(a, 8) for different values of shape parameter a. 


variables. Granados and Sanabria*! studied independence neutrosophic random variables. Neutrosophic has 
many applications in many fields such as decision making, machine learning, intelligent disease diagnosis, com- 
munication services, pattern recognition, social network analysis and e-learning systems, physics, sequences 
spaces and so on for more details see’”*?->! and*. This paper is organized as follows, in section “Neutrosophic 
Birnbaum-Saunders distribution and its statistical properties’, we introduce the new distribution NBS(ay, By). 
and derived its statistical properties. In section “Parameter estimation’, Bayesian and non-Bayesian estimation 
methods are discussed to estimate the parameters of new distribution. In section “Simulation and comparative 
study’, the Monte-Carlo simulation and comparative study is performed to investigate the behavior of different 
estimates for the parameters of our distribution and compare between different estimates of parameters of new 
distribution. In section “Comparative study using real application’, real life data analysis is introduced. Finally, 
in section “Conclusion”, the conclusion of our study is introduced. 


Neutrosophic Birnbaum-Saunders distribution NBS(ay, By). and its statistical 
properties 

In this section, we introduce the new distribution which called neutrosophic Birnbaum-Saunders distribution 
and denoted by NBS(ay, Bn). where @ is a shape parameter and f is the scale parameter. We use Mathematica 
13.1 in all calculations in this section, for more details see*’. 


Definition 2: (The neutrosophic probability density function and neutrosophic cumulative distribution func- 
tion of NBS(an, Bn)) 

Let Iy € (Iz, Iv) be an indeterminacy interval, where N is the neutrosophic statistical number and let 
Xy = X, + XuI be a random variable following neutrosophic Birnbaum-Saunders with scale parameter By 
and shape parameter ay. If the neutrosophic probability density function (npdf) and neutrosophic cumulative 
distribution function (ncdf) are defined as follows respectively, 


exp(=2yss—U" ) By + Bnxn) 


ah) BNXNi 


fn) = (1+ In), xn > 0,aNn > 0, Bn > 0. (4) 
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Note that, the neutrosophic distribution go to the classical distribution when Iy = 0. Figures 4 and 5 show 
npdf and ncdf for different values of wy and By. Two show the changes in the distribution curve when the shape 
parameter takes different values. Also we can see the effect of indeterminacy parameter on curves. 


Definition 3: (The neutrosophic reliability function and neutrosophic hazard function of NBS(an, Bn)) 
The neutrosophic reliability function of Xy is a random variable following neutrosophic Birnbaum-Saunders 
with scale parameter By and shape parameter oy is defined as, 


R(xn) = 1 — F(xy), 
—1+ Byxn 
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1— —({ 1+ Erf} ——_—_ 
5( (|e 


and the neutrosophic hazard function of Xy is defined as, 


(6) 


| )a +n. xn > 0,aNn > 0, By > 0. 
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(7) 
and denoted by nhf. Figure 6 shows the nhf for different values of ay and By. The Figure show the changes in the 
distribution curve when the shape parameter takes different values. Also we can see the effect of indeterminacy 
parameter on curves. 


Now,we discuss some statistical properties of new proposed distribution NBS(an, Bn) such as mode, median, 
moments, moment generating function, quantile function, order statistics, entropy. 

I. Mode: 

To find the mode of neutrosophic Birnbaum- Saunders distribution solve the following nonlinear equation 
with respect to xy, 


exp (=Byau—U") (1 + In)xn (-1 + Bnxn(-1 + 30%, + Bnxn (1 + ay, + Bnxn))) 


ah) BNXNi 
= 0. 
aViva} (Busy) 
Then, the mode at xy = Root[—1 + (—By + 3an7B)#1 + (By? + an? By7)#1? + By? #13 &, 1] 
Where, 0 < ay < Root[—64 + 64#1? — 92#14 + 9#1°&, 2, 0]&&By > 0)||(Root[—64 + 64#1? — 92#14 + 9#1°&, 
1,0] <an < 0&& Bn > 0)||Bn < 0. 
When ay = 1.5, 8 = 2 then xy = 0.0794. 
II. Median 
The median of NBS(an, Bn) is given by. 
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Figure 6. The nhf for different values of ay and By. 
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When ay = 1.5, By = 2,In = 0.2 thenm = (0.3651, 0.6846). 
III. r-th moments of origin 
The r-th moments of origin of NBS(ay, By) is defined as 
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IV. Mean 
The mean of NBS(ay, Bn) is given by, 


A= JF anflaw den, 
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V. Variance 
The Variance of NBS(ay, Bn) is given by. 
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VI. Moment generating function 
The moment generating function of NBS(ay, By) is given by, 
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VII. Characteristic function 
The characteristic function of NBS(ay, By) is given by, 
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where i = /—1. 
VIII. Cumulant generating function 
The cumulant generating function of NBS(ay, By) is given by. 


Cxy (t) = Log[®xy (t)], 
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wherei = /—1. 
IX. Quantile function 
The quantile function of NBS(ay, By) is given by. 


F(Qxy (p)) =P, 
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where, k? = InverseErf [ae] : 
X. Order statistics 
For given any random variables Xn .. . Xwn, the order statistics. Xy1) ... XN) are also random variables, 
defined by sorting the values of Xj . . . Xyw in increasing order. For a random sample Xnv1) .. .Xwqy) the npdf 
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XI. Entropy 


Entropy is considered one of the most popular measures of uncertainty.** introduced the differential entropy 
H(X)as follows: 


H(X) =— / f (x)Log (f (x)) dx 


Rényi> introduced Renyi entropy which finds its source in the information theory. He defined the Renyi 


entropy as follow: 
= 3 
Log / (f(x)) ax : 
6 oo 
Where, 6 4 landé > 0. 
Tsallis°® introduced q-entropy which comes from statistical physics. He defined the q-entropy as follows: 


H3(X) = ylt-F f F (roo). 


Where, q # land q > 0. Now, the three entropies are defined for NBS(ay, Bn) as follows: 
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Parameter estimation 
In this section, maximum likelihood and Bayesian estimation methods were used to estimate the parameters of 
our new proposed distribution neutrosophic Birnbaum-Saunders distribution. 


where, k = 


Maximum likelihood estimation method 
Let Xy1...XNn be a random sample from NBS(ay, Bn). Then the likelihood function is given by. 
n 


L(xyi3 ay, BN) = [] fom: an, Bn). 


i=1 


n 2 
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and its corresponding log-likelihood function is given by, 


1 1 
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Now get the derivatives of log-likelihood function with respect to ay and fy to get the maximum likelihood 
estimates for ay and By which denoted by @y and By as follows: 


dlogLenian, By) _ om | 1 3 (—1+ Byxni)” 
dan ay  an* By XNi 
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Normal equations can’t solve analytic. So, we can't get the closed form for @y and By. Hence, numerical 
method is used to solve these equations. Since the maximum likelihood estimates for unknown parameters of 
new proposed distribution @y and By can’t get in closed form, so the exact distributions of these parameters not 
derived, so we derive the asymptotic confidence intervals of these parameters. For large sample and ay > 0 and 
By > 0. The ay and B Ny are bivariate normal distribution with the mean ay and Ay and covariance matrix I r 
Where I! is the inverse of information matrix, where, 


ma ( Var(an) ae) 
. Cov(@y, Bn) Var(Bn) J” 


For more details see*’. Now the 100(1 — y)% confidence interval of parameters ay and By are 


ant Zy/2W Var(an) and Bw + Zy/2\/ Var (Bn) respectively. 


Bayesian estimation method 

The Bayes estimation using MCMC technique is used to estimate the unknown parameters of new proposed 
distribution NBS(an, Bw). For more details of MCMC technique using Gibb sampling procedure see**”’. Also, 
for more details of MCMC technique using Metropolis Hasting (MH) method see® and°!. The two methods are 
used to generate samples from the posterior density function to compute point Bayes estimators for unknown 
parameters and construct credible intervals. For this aim, we suppose that independent gamma prior distribu- 
tions for unknown parameters of new proposed distribution NBS(an, By) as follows: 
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In this case the joint prior distribution of ay and By is given by, 
(an, BN) = 11 (On )772 (BN), 


And the joint posterior is given by, 
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Under square error loss the Bayes estimates of ay and By are given by. 


m* (an, Bn|Xn ) = 
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These estimates can’t be computed analytically. So, we use MCMC method using MH technique to get the 
aB and B¥ as follows: 


i. Choose initial values of a”) and B i a 


ii. Suppose the values of wy and By at the k" step by a? and pe. 


iii. Generate a) using 2* (co An? Xv ) and z* (Bx leh? Xy ) respectively. 


iv. Repeat step 3 N-times. 
v. Compute Bayes estimates of wy and By as follows: 


where B is the burn-in period. 
vi Compute (100 — y)% HPD credible intervals for ay and By as follows: 


(concyyee(a—z)) (8x4). Bucy): 


Note that we use R-Studio Software to get results in this section using many packages such that, nlme, MASS, 
coda, mcmce, distr, VGAM and RCPP. 


Simulation and comparative study 

In this section, we perform Monte-Carlo simulation study to investigate the behavior of two different estima- 
tors for parameters of new proposed distribution NBS(ay, By) maximum likelihood estimators and Bayesian 
estimates according to different sample sizes, different start values of ay and By and different indeterminacy 
measure. Also, we introduce comparative study to compare between maximum likelihood estimates (MLE’s) 
and Bayesian estimates to get the best for parameters of new proposed distribution NBS(ay, Bn). Compare 
between classical and neutrosophic version of BS distribution to show the flexibility of neutrosophic version. 
Finally, compare Bayesian estimates for different prior distributions. For the aim of comparative study, we use 
bias and mean square error (MSE) to compare between different point estimators. Use also the Akaike informa- 
tion criterion (AIC) to compare maximum likelihood estimator for classical and neutrosophic version but use 
asymptotic confidence length (ACL) to compare between different interval estimators. Now we perform these 
studies according to the following steps: 


i. Choose the different initial values of (wy, By) = (1.25, 3), (0.5, 3), (1, 3). 
ii For Bayesian estimators choose different values of the parameter for gamma prior follows 
(a,b) = (1,2), Ui, 1), 1). 
iii, Use two indeterminacy measure Iy = (0.2, 0.5), (0.6, 0.8). Note that when Iy = 0, the classical version of 
Birnbaum Saunders is obtained. 
iv. Generate different sample sizes n = 50, 100, 200, 500. 
v. Find point estimators using maximum likelihood estimation method and Bayesian estimation method. 
vi. Calculate asymptotic confidence interval and credible interval. 

vii. Perform the comparative study by calculating bias, MSE and AIC to compare MLE’s in both cases classical 
and neutrosophic version of BS distribution, calculating bias, MSE for comparing MLE’s and Bayesian 
estimates in both cases classical and neutrosophic version of BS distribution. And use ACL for comparing 
interval estimation. 


All calculations in this section we use R package For more details about R-Package see®. Results of simulation 
and comparative study between our new proposed distribution neutrosophic Birnbaum-Saunders distribu- 
tion and its classical version Birnbaum-Saunders distribution are shown in Tables 1, 2, 3, 4, 5 and 6 which get 
using R-Studio Software. From These Tables, we get, Tables 1 and 2 shown bias’s, MSE’s and AIC for MLE’s for 
BS(a,B) and NBS(qn, Bn) respectively, we get NBS(an, Bn) has smaller bias’s, MSE’s and AIC than BS(a, B) so, 
we can decided that the neutrosophic version is better than the classical version. Also, in the context of interval 
estimation the NBS(ay, Bn) has smaller ACL for asymptotic confidence interval than BS(a, B). So, For MLE’s 
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BS(a, B) 

a B a B 
n Bias MSE Bias MSE AIC ACL 
(a, B) = (1.25.5, 3) 
50 0.0193 | 0.0188 0.03350 | 0.0563 9963.393 | 0.1654 | 0.5138 
100 0.0088 | 0.0077 0.0165 0.0273 | 32,422.27 | 0.1563 | 0.3786 
200 0.0045 | 0.0041 0.0084 0.0141 | 132,107.5 | 0.0998 | 0.2654 
500 0.0018 | 0.0016 0.0034 0.0060 | 826,131.3 | 0.0593 | 0.1550 
(a, B) = (0.5, 3) 
50 0.0039 | 0.0007 0.0508 0.1292 8980.385 | 0.0744 | 0.1795 
100 0.0018 | 0.0003 0.0239 0.0572 | 34,746.18 | 0.0716 | 0.0167 
200 0.0007 | 0.0001 0.0125 0.0314 | 128,830.6 | 0.0453 | 0.0960 
500 0.0003 | 5.16*10~5 | 0.0050 0.0125 | 810,005.7 | 0.0283 | 0.0618 
(a, B) = (1,3) 
50 0.0138 | 0.0095 0.0428 0.0916 8744.334| 0.1340 | 0.3390 
100 0.0069 | 0.0048 0.0204 0.0419 | 36,114.63 | 0.1169 | 0.2579 
200 0.0034 | 00,023 0.0102 0.0208 | 139,516.5 | 0.0743 | 0.1876 
500 0.0013 | 0.0008 0.0040 0.0080 | 812,134.2 | 0.0516 | 0.1246 


Table 1. Bias’s, MSE’s, AIC’s and ACLs for MLE’s for BS(q, B). 


the NBS(an; Bn) has good behavior than BS(a, 8). Also, we get for two versions bias’s and MSE’s decrease when 
sample size increase. For Bayesian estimator results are shown in Tables 3, 4, 5 and 6, we get also, NBS(an; Bn) 
has smaller bias’s and MSE’s than BS(q, 8) as shown in Tables 5 and 6. Also, for credible intervals NBS(an, Bn) 
has smaller ACL than BS(q@, 8). The behavior of Bayesian estimation got for different three prior distributions. 
So, we can decide also, NBS(a; Bn) has good behavior than BS(a, B). 


Comparative study using real application 
The main aim of this section is a comparative study between the NBS(ay, Bn) and BS(q, 6). We introduced two 
real applications as follows: 


Application 1 
Based on data of alloy melting points. For more details see® which mentioned that A combination of material 
constituents, including at least one metal, makes up an alloy. In general, evaluating melting points is quite chal- 
lenging, therefore observations are indeterministic and can be communicated in intervals as follows: 

[563.3, 545.5], [529.4, 511.6], [523.1, 503.5], [470.1,449.2], [506.7, 489.0], [495.6, 479.1], [495.3, 467.9],[520.9, 
495.6], [496.9, 472.8], [542.9, 519.1], [505.4,484.0], [550.7, 525.9], [517.7, 500.9], [499.2, 483.0],[500.6, 480.0], 
[516.8, 499.6], [535.0, 515.1], [489.3,464.4]. 


Application 2 
The data represent the lifetime of batteries. The lifetime in 100hours of 23 batteries is given as: 

[2.9,3.99], [5.24,7.2],[6.56,9.02], [7.14,9.82], [11.6,15.96], [12.14,16.69], [12.65,17.4], [13.24,18.21], 
[13.67,18.79], [13.88,19.09], [15.64,21.51], [17.05,23.45], [17.4,23.93], [17.8,24.48], [19.01,26.14], [19.34,26.59], 
[23.13,31.81], [23.34,32.09],[26.07,35.84], [30.29,41.65], [43.97,60.46], [48.09,66.13], [73.48,98.04]. For more 
details see®* 

Figure 7 shows the architectural diagram of the proposed algorithm in this section. To show the performance 
of our new distribution the NBS(ay,, By). We compare it with its classical version BS(q, 6) using three statistical 
criteria log -likelihood (— 2LL), Akaike’s Information Criteria (AIC) and Bayesian Information Criteria (BIC), 
where, 


AIC = —2LogL(@) + 2k, BIC = —2LogL(©) + kLog{[n]. 


where, ©: the vector of distribution parameters, L| © |}: the likelihood function, k: the number of estimates, n: 


the data size. The small value of -2LL, AIC and BIC mean good-fit distribution. Table 7 shows the result of com- 
parison between our new distribution and other distribution under classical statistics. The result in Table 7 shows 
our new distribution is better for this data than its classical version in two applications because all goodness of 
fit tests having smaller values in neutrosophic version than the classical version. 


Conclusion 

A new distribution introduced which called neutrosophic Birnbaum-Saunders distribution and denoted by 
NBS(an, Bn). some statistical properties such as neutrosophic probability density function, neutrosophic cumu- 
lative distribution function, neutrosophic hazard function, neutrosophic mean, mode, median, variance, moment 
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n 


NBS(an, Bn) 
ON Bn oN Bn 
Bias MSE Bias MSE AIC ACL 


(an; Bn) = (1.25, 3), Iy = (0.2, 0.5) 


50 (0.0093, 0.0153) (0.0044, 0.0118) (0.0235, 0.0295) | (0.0278, 0.0437) —_| (9922.846,9945.161) | (0.1651,0.1652) —_| (0.5118, 0.5128) 
100 (0.0038, 0.0068) (0.0014, 0.0046) (0.0115, 0.0145) | (0.0133, 0.0211) _| (32,341.18, 32,385.81) | (0.1556, 0.1562) —_ | (0.3754, 0.3767) 
200 (0.0020, 0.0035) (0.0008, 0.0025) (0.0059, 0.0074) | (0.0070, 0.0110) —_| (131,945.3, 132,034.5) | (0.0005, 0.0007) | (0.2492, 0.2653) 
500 (0.0008, 0.0014) (0.0003, 0.0010) (0.0024, 0.0030) | (0.0030, 0.0047) —_| (825,725.8, 825,949) —_| (0.0587, 0.0590) _| (0.14831, 0.1524) 


(an; Bn) = (0.5, 3), In = (0.2, 0.5) 


50 (- 0.0060 , — 2.81075) | (4.09*10-8, 0.0006 ) _| (0.0408, 0.0468) —_| (0.0834, 0.1097) _| (8939.838, 8962.153) | (0.0740, 0.0743) __| (0.1782, 0.1792) 
100 (— 0.0003, 0.0001) —_| (2.6*10~*, 0.0002) (0.0189, 0.0219) | (0.0358, 0.0481) —_| (34,665.09, 34,709.72) | (0.0710, 0.0713) | (0.0162, 0.0163) 
200 (- 0.0002, 0.0017) | (9.01*10~, 0.0001) (0.0100, 0.0115) | (0.0201, 0.0266) _| (128,668.4, 128,757.6) | (0.0450, 0.0452) —_| (0.0952, 0.0958) 
500 (- 7.8*10-°, — 0.0006) | (3.08*10~°, 0.0001) (0.0040, 0.0046) | (0.0080, 0.0106) —_| (809,600.2, 809,823.4) | (0.0280, 0.0282) | (0.0613, 0.0614) 


(an, Bn) = (1,3), In = (0.2,0.5) 


50 (0.0038, 0.0098) (0.0007, 0.0048) (0.0328, 0.0388) | (0.0538, 0.07534) | (8703.788, 8726.102) | (0.1338, 0.1339) _| (0.3367, 0.3378) 
100 (0.0019, 0.0049) (0.0003, 0.0024) (0.0154, 0.0184) | (0.0239, 0.0341) —_| (36,033.54, 36,078.17) | (0.1156, 0.01159) _ | (0.2558, 0.2569) 
200 (0.0009, 0.0024) (0.0001, 0.0011) (0.0077, 0.0062) | (0.0118, 0.0169) —_| (139,354.3, 139,443.5) | (0.0740, 0.0742) —_| (0.1816, 0.1845) 
500 (0.0009, 0.0003) (5.45*10—°, 0.0004) (0.0030, 0.0036) | (0.0045, 0.0065) —_| (811,728.7, 811,951.8) | (0.0512,0.0514) | (0.1241, 0.1245) 


(an, Bn) = (1.25, 3), Iv = (0.6, 0.8) 


50 (— 0.0001, 0.0038) (8.32*1079, 0.0007) (0.0154, 0.0190) —_| (0.0120, 0.0190) _| (5667.594, 5679.373) | (0.0819, 0.0823) _| (0.4594, 0.4646) 
100 (0.0009, 0.0029) (8.46*1075, 0.0008) (0.0101, 0.0121) | (0.0103, 0.0148) —_| (30,906.83, 30,930.38) | (0.0395, 0.0396) | (0.3249, 0.3290) 
200 (0.0005, 0.0006) (5.15*107°, 0.0004) (0.0044, 0.0054) | (0.0040, 0.0060) —_| (122,478.9, 122,526) | (0.0017, 0.0022) —_| (0.2528, 0.2529) 
500 (0.0002, 0.0006) (2.84*10-°, 0.0002) (0.0018, 0.0022) | (0.0017, 0.0025) —_| (844,985.4, 845,573.2) | (0.0401, 0.0410) | (0.1425, 0.1513) 


(an, Bn) = (0.5, 3), In = (0.6, 0.8) 


50 (- 1.0131, - 0.0091) —_| (0.0041, 0.0086) (0.0341, 0.0381) | (0.0584, 0.0729) —_| (7632.389, 7644.168) | (0.0640, 0.0641) _| (0.1620, 0.1622) 
100 (— 0.0063, — 0.0043) —_| (0.0019, 0.0040) (0.0168, 0.0188) | (0.0284, 0.0355) —_| (32,526.3, 32,549.85) | (0.0656, 0.0658) _| (0.1423, 0.1426) 
200 (— 0.0031, - 0.0021) —_| (0.0009, 0.0019) (0.0082, 0.0092) | (0.0137,0.0172) —_| (134,481.5, 134,528.6) | (0.0422, 0.0433) | (0.0051, 0.0054) 
500 (- 0.0012, - 0.0008) —_| (0.0003, 0.0007) (0.0033, 0.0037) | (0.0057,0.0071) —_| (840,294.8, 840,412.6) | (0.0281, 0.0282) —_| (0.0523, 0.0532) 


(ay, Bx) = (1,3), In = (0.6, 0.8) 


50 (— 0.0025, 0.0014) (0.0001, 0.0033) (0.0255, 0.0295) | (0.0325, 0.0435) —_| (8349.03, 8360.808) —_| (0.1208, 0.1209) _| (0.3618, 0.3647) 
100 (— 0.0019, 6.38*10->) | (4.07*10~7, 0.0003) (0.0126, 0.0146) | (0.0159,0.0214) —_| (27,309.47, 27,333.03) | (0.1107, 0.1109) | (0.2479, 0.2499) 
200 (— 0.0007, 0.0002) (1.60*10-5, 0.0001) (0.0060, 0.0070) | (0.0073, 0.0099 (129,253.6, 129,300.7) | (0.0630, 0.0632) | (0.1808, 0.1824) 
500 (— 0.0002, 0.0001) (6.08*10~®, 4.19*10-5) | (0.0024, 0.0028) —_| (0.0030, 0.0041) —_| (800,043, 800,160.8) | (0.0513, 0.0515) —_‘| (0.1178, 0.1244) 


Tablae 2. Bias’s, MSE’s, AIC’s and ACLs for MLE’s for NBS(q, f). 


Credible intervals 

BS(a, B) 

(a, B) = (1.25, 3) 

ACL 

Prior 1 Prior 2 Prior3 
n a B a B a B 
50 0.9457 | 0.9439 | 0.9533 | 0.9322 | 0.9385 | 0.9589 
100 0.9411 | 0.9545 | 0.9462 | 0.9315 | 0.9471 | 0.9604 
200 0.9560 | 0.9337 | 0.9608 | 0.9458 | 0.9594 | 0.9474 
500 0.9522 | 0.9245 | 0.9362 | 0.9513 | 0.9549 | 0.9313 
(a, B) = (0.5, 3) 
50 0.9406 | 0.9413 | 0.9480 | 0.9517 | 0.9443 | 0.9607 
100 0.9448 | 0.9422 | 0.9296 | 0.9505 | 0.9416 | 0.9395 
200 0.9429 | 0.9313 | 0.9281 | 0.9530 | 0.9251 | 0.9441 
500 0.9466 | 0.9486 | 0.9489 | 0.9465 | 0.9353 | 0.9521 
(a, B) = (1,3) 
50 0.9279 | 0.9498 | 0.9489 | 0.9378 | 0.9442 | 0.9456 
100 0.9531 | 0.9385 | 0.9372 | 0.9425 | 0.9439 | 0.9499 
200 0.9330 | 0.9340 | 0.9389 | 0.9498 | 0.9382 | 0.9372 
500 0.9423 | 0.9378 | 0.9411 | 0.9414 | 0.9513 | 0.9437 


Table 3. The ALC’s for Credible intervals for BS(a, B). 
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Credible intervals 


NBS(q, B) 


(on, Bn) = (1.25, 3), Iv = (0.2, 0.5) 


ACL 


Prior 1 


Prior 2 


Prior3 


oN 


By 


oN 


By 


aN 


By 


(0.9299, 0.9450) 


(0.9406, 0.9420) 


(0.9321, 0.9329) 


0.9276, 0.9238) 


(0.9286, 0.9295) 


(0.9466, 0.9538) 


(0.9389, 0.9391) 


(0.9380, 0.9466) 


(0.9435, 0.9446) 


0.9257, 0.9280) 


(0.9465, 0.9473) 


(0.9491, 0.9578) 


(0.9362, 0.9387) 


(0.9227, 0.9324) 


(0.9513, 0.9542) 


0.9389, 0.9401) 


(0.9584, 0.9586) 


(0.9213, 0.9456) 


500 


(0.9402, 0.9419) 


(0.9212, 0.9222) 


(0.9391, 0.9579) 


( 
( 
( 
( 


0.9271, 0.9403) 


(0.9464, 0.9540) 


(0.9240, 0.9262) 


(an, Bn) = (0.5, 3), Iy = (0.2, 0.5) 


50 


0.9322, 0.9518) 


(0.9386, 0.9393) 


(0.9433, 0.9478) 


0.9433, 0.9487) 


(0.9400, 0.9441) 


(0.9416, 0.9492) 


100 


0.9296, 0.9317) 


(0.9338, 0.9344) 


(0.9246, 0.9267) 


0.9339, 0.9468) 


(0.9362, 0.9415) 


(0.9353, 0.9392) 


200 


0.9192 0.9337) 


(0.9228, 0.9234) 


(0.9195, 0.9267) 


0.9446, 0.9461) 


(0.9194, 0.9212) 


(0.9402, 0.9520) 


500 


0.9402, 0.9417) 


(0.9333, 0.9458) 


(0.9333, 0.9437) 


( 
( 
( 
( 


0.9378, 0.9451) 


(0.9307, 0.9341) 


(0.9426, 0.9433) 


(an, By) = (1,3), In = (0.2, 0.5) 


50 


(0.9146, 0.9203) 


(0.9354, 0.9362) 


(0.9471, 0.9485) 


(0.9382, 0.9433, 


(0.9323, .9335) 


(0.9348, 0.9500) 


100 


(0.9432, 0.9434) 


(0.9339 0.9384) 


(0.9363, 0.9364) 


(0.9315, 0.9436 


(0.9247, 0.9406) 


(0.9371, .9390) 


200 


(0.9287, 0.9318) 


(0.9293, 0.9301) 


(0.9337, 0.9349) 


(0.9377, 0.9404 


(0.9125, 0.9232) 


(0.9362, 0.9349) 


500 


(0.9349, 0.9394) 


(0.9312, 0.9374) 


(0.9262, 0.9320) 


) 
) 
) 
) 


(0.9430, 0.9451 


(0.9313, 0.9340) 


(0.9304, 0.9366) 


(ay, Bn) = (1.25, 3), In = 


(0.6, 0.8) 


50 


(0.9443 0.9446) 


(0.9301, 0.9409) 


(0.9355, 0.9488) 


(0.9285, 0.9293) 


(0.9250, 0.9295) 


(0.9457, 0.9576) 


100 


(0.9250, 0.9316) 


(0.9391, 0.9446) 


(0.9356, 0.9423) 


(0.9289, 0.9298) 


(0.9407, 0.9462) 


(0.9351, 0.9597) 


200 


(0.9351, 0.9434) 


(0.9265, 0.9274) 


(0.9368, 0.9492) 


(0.9443, 0.9457) 


(0.9270, 0.9319) 


(0.9314, 0.9409) 


500 


(0.9382, 0.9439) 


(0.9175, 0.9201) 


(0.9325, 0.9344) 


(0.9458, 0.9462) 


(0.9418, 0.9434) 


(0.9232, 0.9274) 


(ay, Bn) = (0.5, 3), Iv = (0.6, 0. 


8) 


50 


(0.9314, 0.9323) 


(0.9335, 0.9373) 


(0.9249, 0.9380) 


0.9458, 0.9461) 


(0.9370, 0.9378) 


(0.9369, 0.9605) 


100 


(0.9342, 0.9444) 


(0.9416, 0.9420) 


(0.9137, 0.9215) 


0.9470, 0.9482) 


(0.9360, 0.9377) 


(0.9344, 0.9394) 


200 


(0.9309, 0.9318) 


(0.9300, 0.9302) 


(0.9162, 0.9204) 


0.9442, 0.9448) 


(0.9163, 0.9222) 


(0.9416, 0.9433) 


500 


(0.9455, 0.9464) 


(0.9471, 0.9484) 


(2.9324, 2.9396) 


( 
( 
( 
( 


0.9245, 0.9375) 


(0.9348, 0.9350) 


(0.9394, 0.9443) 


(ay, By) = (1, 3), In = (0.6, 0.8) 


50 


(0.9220, 0.9259) 


(0.9456, 0.9464) 


(0.9318, 0.9420) 


0.9238, 0.9366) 


(0.9370, 0.9389) 


(0.9381, 0.9433) 


100 


(0.9523, 0.9524) 


(0.9323, 0.9341) 


(0.9178, 0.9249) 


0.9338, 0.9411) 


(0.9287, 0.9410) 


(0.9452, 0.9460) 


200 


(0.9311, 0.9321) 


(0.9314, 0.9321) 


(0.9334, 0.9373) 


0.9373, 0.9446) 


(0.9272, 0.9362) 


(0.9328, 0.9352) 


500 


(0.9383, 0.9420) 


(0.9356, 0.9366) 


(0.9335, 0.9347) 


( 
( 
( 
( 


0.9305, 0.9405) 


(0.9381, 0.9465) 


(0.9403, 0.9429) 


Table 4. The ALC’s for credible intervals for NBS(q, B). 


BS(a, B) 
Prior 1 Prior 2 Prior3 
a B a B a B 
n Bias MSE Bias MSE Bias MSE Bias MSE Bias MSE Bias MSE 
(a, B) = (1.25, 3) 
50 0.0119 | 0.0075 | 0.0175 | 0.0306 | 0.0144 | 0.0131 | 0.0172 | 0.0316 | 0.0076 | 0.0062 | 0.0135 | 0.0106 
100 | 0.0096 | 0.0446 | 0.0087 | 0.0153 | 0.0080 | 0.0314 | 0.0085 | 0.0163 | 0.0036 | 0.0013 | 0.0088 | 0.0123 
200 | 0.0032 | 0.0030 | 0.0043 | 0.0076 | 0.0042 | 0.0209 | 0.0041 | 0.0066 | 0.0015 | 0.0010 | 0.0013 | 0.0036 
500 | 0.0058 | 0.0366 | 0.0017 | 0.0030 | 0.0014 | 0.0011 | 0.0015 | 0.0027 | 0.0021 | 0.0022 | 0.0014 | 0.0003 
(a, B) = (0.5, 3) 
50 0.0084 | 0.0062 | 0.0175 | 0.0306 |-0.0008 | 0.0018 | 0.0250 | 0.0625 | -0.0040 | 0.0064 | 0.0030 | 0.0070 
100 | 0.0086 | 0.0221 | 0.0087 | 0.0153 | 0.0035 | 0.0013 | 0.0125 | 0.0312 |-0.0050 | 0.0100 | 0.0050 | 0.0212 
200 | 0.0041 | 0.0080 | 0.0043 | 0.0076 | 0.0011 | 0.0005 | 0.0062 | 0.0156 | 7.8%10~7| 0.0002 | 0.0033 | 0.0125 
500 0.0025 | 0.0081 | 0.0017 | 0.0030 |-0.0011 | 0.0008 | 0.0025 | 0.0006 | — 0.0013 0.0037 | 0.0012 | 0.0003 
(a, B) = (1,3) 
50 0.0135 | 0.0109 | 0.0200 | 0.0400 0.0046 | 0.0094 | 0.0019 | 0.0038 0.0027 0.0073 | 0.0150 | 0.0425 
100 | 0.0030 | 0.0007 | 0.0100 | 0.0200 | 0.0028 | 0.0010 | 0.0009 | 0.0019 | 0.0051 | 0.0049 | 0.0075 | 0.0212 
200 | 0.0024 | 0.0012 | 0.0050 | 0.0100 | 0.0041 | 0.0034 | 0.0048 | 0.0018 |-0.0007 | 0.0002 | 0.0038 | 0.0106 
500 0.0016 | 0.0014 | 0.0020 | 0.0040 0.0009 | 0.0006 | 0.0018 | 0.0038 | — 0.0041 0.0187 | 0.0015 | 0.0043 
Table 5. Bias’s and MSE’s for Bayesian estimation for BS(q, f). 
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NBS(q, B) 
Prior 1 Prior 2 Prior3 
oN Bu oN Bu ON By 
n Bias MSE Bias MSE Bias MSE Bias MSE Bias MSE Bias MSE 
(ay, BN) = (1.25, 3), In = (0.2, 0.5) 
- (0.0064, (0.0023, (0.0075 (0.0180, (0.0041, (0.0018, (0.0065, (0.0181, (— 0.0069, — | (0.0034, (0.0126, (0.0006, 
0.0088) 0.0039) 0.0115) 0.0244) 0.0134) 0.0129) 0.0135) 0.0243) 0.0063) 0.0100) 0.0130) 0.0009) 
100 _ | (0.0014; (0.0011, (0.0038, (0.0091, (0.0037, (0.0014, (0.0037, (0.0090, (0.0091, (0.0007, (0.0039, (0.0113, 
0.0021) 0.0031) 0.0067) 0.0122) 0.0047) 0.0022) 0.0068) 0.0122) 0.0015) 0.0012) 0.0086) 0.0121) 
209 | (0.0169, — | (0.0024, (0.0019. (0.0045, (— 0.0016, | (0.0006, (0.0018, (0.0043, (0.0011, (0.0004, (0.0002, (0.0004, 
0.0028) 0.0029) 0.0031) 0.0061) 0.0016) 0.0052) 0.0034) 0.0062) 0.0022) 0.0009) 0.0009) 0.0006) 
500 _ | (0.0001, (0.0005, (0.0008, (0.0018, (0.008, (0.0003, (0.0007, (0.0019, (- 0.0007, | (1.6*10~5 (0.0012, (0.0001, 
0.0006) 0.0017) 0.0014) 0.0024) 0.0009) 0.0004) 00,013) 0.0025) 0.0020) 0.0012) 0.0013) 0.0003) 
(ay, Bn) = (0.5, 3), In = (0.2, 0.5) 
a (— 0.0334 — | (0.0003, (0.0150, (0.0300, (— 0.0356, — | (0.0009, (0.0020, (0.0522 (— 0.0139, — | (0.0091, (0.0025, (0.0062, 
0.0016) 0.0004) 0.0170) 0.0303) 0.0084) 0.0015) 0.0250) 0.0625) 0.0133) 0.0043) 0.0330) 0.0067) 
100 | (9.0093, - | (0.0064, (0.0075, (0.0113, (— 0.0103, — | (0.0010, (0.0115, (0.0212, (— 0.0177, — | (0.0006, (0.0025, (0.0112, 
0.0063) 0.0096) 0.0085) 0.0127) 0.0102) 0.0012) 0.0117) 0.0302) 0.0004) 0.0009) 0.0035) 0.0147) 
209 | ©.0-0099, - | (0.0040, (0.0038, (0.0066, (— 0.0046, — | (0.0001, (0.0053, (0.0106, (- 0.0114, | (0.0001, (0.0013, (0.0100, 
0.0073) 0.0074) 0.0042) 0.0073) 0.0101) 0.0004) 0.0057) 0.0146) 0.0013) 0.0002) 0.0025) 0.0114) 
500 | 90-0023, - | (0.0008, (0.0012, (0.0025, (— 0.0011, — | (0.0002, (0.0015, (0.0002, (— 0.0028, — | (0.0033, (0.0006, (0.0001, 
0.0011) 0.0027) 0.0016) 0.0030) 0.0005) 0.0007) 0.0023) 0.0006) 0.0021) 0.0036) 0.0010) 0.0002) 
(ay; By) = (1,3), Iy = (0.2, 0.5) 
50 (9.40*10—5, | (0.0024, (0.0100, (0.0250, (0.0023, (0.0004, (0.0019, (0.0009, (— 0.0418, — | (0.0022, (0.0001, (0.0200, 
0.0081) 0.0054) 0.0160) 0.0328) 0.0039) 0.0018) 0.0040) 0.0018) 0.0329) 0.0026) 0.00130) 0.0400) 
100 | (9.0048, - | (0.0001, (0.0050, (0.0125, (— 0.0038, — | (0.0015, (0.0009, (0.0002, (— 0.0020, | (0.0009, (0.0008, (0.0100, 
0.0026) 0.0007) 0.0080) 0.0164) 0.0036) 0.0021) 0.0010) 0.0005) 0.0004) 0.0013) 0.0060) 0.0200) 
200 | (0.0042, - | (0.0006, (0.0025, (0.0065, (5.49*10—5, | (0.0002, (0.0025, (0.0030, (— 0.0034, — | (0.0001, (0.0004, (0.0050, 
0.0002) 0.0011) 0.0040) 0.0082) 0.0001) 0.0031) 0.0030) 0.0035) 0.0100) 0.0002) 0.0020) 0.0100) 
500 | (0.0021, — | (0.0005, (0.0010, (0.0025, (— 0.0026, (5.6*10~°, | (0.0001, (0.0010, (— 0.0023, — | (0.0006, (0.0005, (0.0031, 
0.0004) 0.0012) 0.0016) 0.0033) 4.82*10->) | 0.0099) 0.0002) 0.0015) 0.0008) 0.0034) 0.0019) 0.0040) 
(ay, BN) = (1.25, 3), In = (0.6, 0.8) 
50 (— 0.0116, — | (0.0010, (0.0012, (0.0154, (— 0.0293, — | (0.0060, (0.0145, (0.0306, (— 0.0026, | (0.0041, (0.0125, (0.0096, 
0.0026) 0.0028) 0.0055) 0.0168) 0.0010) 0.0067) 0.0165) 0.0308) 0.0028) 0.0045) 0.0130) 0.0104) 
100 | (9.0236, | (0.0033, (0.0008, (0.0077, (- 0.0171, | (0.0011, (0.0067, (0.0153, (— 0.0051, | (0.0007, (0.0080, (0.0114, 
0.0003) 0.0086) 0.0028) 0.0084) 0.0032) 0.0254) 0.0078) 0.0159) 0.0032) 0.0013) 0.0087) 0.0117) 
200 | (0.0041, | (0.0017, (0.0004, (0.0038, (— 0.0092, — | (0.0113, (0.0040, (0.0056, (— 0.0005, — | (6.2*1075 (0.0007, (0.0027, 
0.0029) 0.0026) 0.0014) 0.0042) 0.0059) 0.0157) 0.0041) 0.0064) 0.0036) 0.0004) 0.0012) 0.0033) 
500 | 0.0014, (4.48*10—5, | (0.0001, (0.0015, (- 0.0021, | (0.0002, (0.0008, (0.0017, (— 0.0031, | (0.0014, (0.0010, (0.0001, 
0.0001) 0.0027) 0.0006) 0.0017) 0.0012) 0.0007) 0.0013) 0.0020) 0.0008) 0.0020) 0.0013) 0.0003) 
(ay, Bx) = (0.5, 3), Iy = (0.6, 0.8) 
eh (— 0.02901, | (0.0056, (0.0090, (0.0253, (— 0.0169, — | (0.0009, (0.0150, (0.0530, (— 0.0218, | (0.0052, (0.0020, (0.0005, 
— 0.0283) 0.0059) 0.0130) 0.0297) 0.0105) 0.0018) 0.0170) 0.0607) 0.0122) 0.0054) 0.0025) 0.0007) 
100 | (0.0093, - | (0.0083, (0.0045, (0.0126, (— 0.0300, — | (0.0005, (0.0110, (0.0213, (— 0.0074, | (0.0039, (0.0021, (0.0112, 
0.0008) 0.0088) 0.0065) 0.0148) 0.0089) 0.0008) 0.0125) 0.0288) 0.0003) 0.0057) 0.0035) 0.0115) 
200 | &.0.0148, - | (0.0076, (0.0023, (0.0058, (— 0.0048, — | (0.0001, (0.0022, (0.0126, (— 0.0127, | (0.0001, (— 0.0017, | (0.0084, 
0.0065) 0.0080) 0.0033) 0.0069) 0.0021) 0.0004) 0.0042) 0.0140) 0.0014) 0.0002) 0.0022) 0.0111) 
500 | .0-0034, - | (0.0027, (0.0001, (0.0025, (— 0.0020, — | (0.0004, (0.0015, (0.0002, (— 0.0013, — | (0.0003, (-4*10-4, | (0.0001, 
0.0021) 0.0067) 0.0013) 0.0029) 0.0009) 0.0007) 0.0025) 0.0004) 0.0008) 0.0009) 0.0011) 0.0003) 
(an, Bx) = (1,3), In = (0.6, 0.8) 
50 (— 0.0262, — | (0.0007, (0.0010, (0.0232, (— 0.0381, — | (0.0011, (0.0014, (0.0030, (— 0.0401, | (0.0054, (— 0.0120, | (0.0272, 
0.0010) 0.0108) 0.0019) 0.0208) 0.0124) 0.0064) 0.0019) 0.0035) 0.0058) 0.0059) 0.0019) 0.0039) 
100 | (9.0063, - | (0.0006, (0.0004, (0.0116, (— 0.0065, | (0.0003, (0.0005, (0.0012, (— 0.0363, —_| (0.0033, (— 0.0120, | (0.0018, 
0.0002) 0.0007) 0.0009) 0.0104) 0.0005) 0.0009) 0.0006) 0.0017) 0.0139) 0.0044) 0.0080) 0.0173) 
200 | ©.9-0058, - | (0.0009, (0.0010, (0.0052, (— 0.0033, — | (0.0019, (0.0030, (0.0010, (— 0.0039, | (0.0001, (— 0.003, (0.0068, 
0.0055) 0.0011) 0.0020) 0.0058) 0.0017) 0.0025) 0.0045) 0.0018) 0.0017) 0.0002) 0.0050) 0.0100) 
500 | (.0-0043, - | (0.0003, (0.0001, (0.0021, (— 0.0032, — | (0.0001, (0.0012, (0.0021, (— 0.0029, | (0.0021, (— 0.0012, | (0.0023, 
0.0003) 0.0013) 0.0004) 0.0023) 0.0007) 0.0004) 0.0018) 0.0025) 0.0002) 0.0044) 0.0020) 0.0040) 


Table 6. Bias’s and MSE’s for Bayesian estimation for NBS(q, B). 


Scientific Reports| (2024) 14:6955 | https://doi.org/10.1038/s41598-024-57438-8 nature portfolio 


www.nature.com/scientificreports/ 


Crisp Input Melting 
points 


Maximum Likelihood Estimation 


Crisp output 


Figure 7. Architectural Diagram of the Proposed Algorithm. 


Application Distribution | MLE - 2LL AIC BIC 


ay = (0.0019, 0.0020) 
By = (0.0044, 0.0047) 


a = 0.0019 


NBS(ay, By) (2929.146,2999.523) | (5860.827,5897.769) | (5864.169,5904.826) 


Application 1 


BS(a, B) ate 11,939,523 23,855.62 23,887.02 
an = (0.0626, 0.0456) 
NBS(ary, By) | BN = (o.78g3; 7546) | (4403:694,4599.646) | (4407.694,4503) (4409.474,4505.427) 
Application 2 
BS(q, ) 7 7 cee 4671.254 4675.254 4677.035 


Table 7. The result of comparison between our new distribution NBS(an, Bn) and BS(q, B). 


of origin, moment generating function, characteristic function, quantile function, cumulant generating function, 
order statistic and entropies. The neutrosophic maximum likelihood estimators and neutrosophic Bayesian esti- 
mators are derived. The simulation study is performed to study the behavior of different estimators at different 
parameter values and different sample sizes. Also, compare between neutrosophic Birnbaum-Saunders distribu- 
tion and its classical version using statistical criteria such as bias, MSE, AIC and ACL. Finally, real data of the 
melting point of alloy and real data of the lifetime of batteries are used to show the validity of NBS(an, By) in real 
life. Also, compare the performance of NBS(ay, By) and BS(a, 6) on this data using AIC, BIC and — 2LL, which 
shows the good performance of NBS(ay, By) than BS(a, 8). The neutrosophic field has more points need more 
and more search so, in future we try to introduce more netrosophic distribution to used solve and describe more 
real life applications also we will try to introduce many researches in different points in neutrosophic statistics. 


Data availability 
The data is given in the paper. 
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